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Abstract 

Exact and asymptotic formulae are displayed for the coefficients A n used in Li's criterion for 
the Riemann Hypothesis. For n — > oo we obtain that if (and only if) the Hypothesis is true, 
A n ~ n(Alogn + B) (with A > and B explicitly given, also for the case of more general zeta 
or L-functions); whereas in the opposite case, \ n has a non-tempered oscillatory form. 

O ! 

Li's criterion for the Riemann Hypothesis (RH) states that the latter is true if and 
only if a specific real sequence {X n } n =i,2,... has all its terms positive jT?l Here we show 
that it actually suffices to probe the A n for their large-n behavior, which fully encodes the 
Riemann Hypothesis by way of a clear-cut and explicit asymptotic alternative. To wit, we 
first represent A n exactly by a finite oscillatory sum Q , then by a derived integral formula 
()12j) . which can finally be evaluated by the saddle-point method in the n — > +oo limit. As 
a result, A n takes one of two sharply distinct and mutually exclusive asymptotic forms: 
if RH is true, A n will grow tamely according to (fTTj) ; if RH is false, A n will oscillate with 
an exponentially growing amplitude, in both + and — directions, as described by (fTHl) . 
This dichotomy thus provides a sharp criterion of a new asymptotic type for the Riemann 
Hypothesis (and for other zeta-type functions as well, replacing (fTTj) by (jT5j) ). 

This work basically reexposes our results of April 2004 announced in |28! , but with an 
uncompressed text; we also update the references and related comments: for instance, we 
now derive as (|24|) a large-n expansion surmised by Maslanka [20] in the meantime. 



*Also at: Institut de Mathematiques de Jussieu-Chevaleret (CNRS UMR 7586), Universite Paris 7, 
F-75251 Paris CEDEX 05, France. 
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1 Background and notations. 

We study the sequence O Ej (in the notations of Li, whose \ n are n times Keiper's) 

A n = [1 - (1 - Vp)1 (n=l,2 l ...), (1) 

p 

where p are the nontrivial zeros of Riemann's C( s ); grouped by pairs within summations 
and products as 

{p — ~ ± iT fc } fc=12v .., Re r fc positive and non-decreasing; (2) 

we also parametrize each such pair by the single number Xk = p(l — p) = \ + t*, 2 . 

We will use the completed zeta function 3(s) (normalized as 2(0) = 2(1) = 1) and a 
symmetrized form of its Hadamard product formula [7| l2"o]. 



x r 8(1-8) 



(3) 



k=i 

we will also use a "secondary" zeta function built over the Riemann zeros, 



Z(a) = J2^ a , Rea>|, (4) 



fc=i 



which extends to a meromorphic function in C having all its poles at the negative half- 
integers, plus one pole at a = +| [T3] of polar part 



Z(\ + e) = ft_ 2 e~ 2 + R_ x a' 1 + O(l) £ _> , (5) 

with 

R_ 2 = (8iry l , iL.i = -(47r) -1 log27r in this case. (6) 



Our results [22] mainly relate to those of Keiper ^2], of which we only learned later 
(thanks to K. Maslanka; they were almost never cited), of Bombieri-Lagarias [2] on Li's 
criterion [T7j , and of Maslanka [2D] • Other earlier works considering the A n are [TJ [3] • 

2 New exact forms for A n . 

To reexpress the A n , we start from their generating function fF2\ IT7] 
Now the infinite product formula in (JHJ) implies 



then, expanding (1 — z) 2j ' by the generalized binomial formula, reordering in powers of 
z and substituting the output into (J7j), we get as first result 



2j 



Other sums related to the Z(k) are ifj = J2 P P * [221 EHl EH| (often denoted Oj, but 
here we use a as variable). It was already known that A n = YS=i(- l ) i+X Q)&i O 
Equation (27)], and that the 3?j in turn are complicated polynomials in the Stieltjes 
constants {~)k}k<j [22] (for X n and 7^ see also [211 EH E] and references therein). Now the 
latter relations boil down to = 1 — (1 — 2~*)C(J) + (— [2EI Equation (46)] simply 
by promoting logarithmic coefficients rjj [UJ|2] (cf. also [TBI Equation (12)]) 



OO , s 

(-S) 



log[<(l + a)] = - E^-V (10) 

n=l 

in place of 

*(! + <> = l-2>~. ■(„_!),■ 

n=l v 7 

The A n thus express as affme combinations of the rjj [21 thm 2]. 
Remarks: 

- the ?•/.,• are the Stieltjes [constants'] cumulants, up to some relabelings [231211; 

- the rjj admit an arithmetic expression over the primes [2, Equation (4.1)]; see also 
[TO] , which cites [23] for the case rj = —7; 

- relation Q can be inverted also in closed form, by the same technique as for [21)1 
Equation (48)]: 
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Z(j) = YX-l) n+1 [ r J ^. (11) 



2j 



J - 71 



n=l 

The expression © for A n has some distinctive advantages: it involves the functional 
equation H(s) = 5(1 — s) through (j3J); and unlike the if}, the Z(j) are positive and gently 
varying factors: the function Z(a) is regular and very smooth for real a > 1. Still, (0) is 
an oscillatory sum, hence difficult to control directly. 

Now an integral representation, equivalent to Q simply by residue calculus, will nev- 
ertheless prove much more flexible: 

(-l) n ni / . , , . , T(a + n)T(a -n) . . 

An = ^ f 1(a) da, 1(a) = 1 r J Z ^ ( 12 ) 

71 J c 1(2CT + 1J 

where C is a positive contour encircling just the subset of poles a = +1, . . . , +n of the 
integrand 1(a). 
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3 Asymptotic alternative for A n , n — > oo. 

The integral formula (|T2*)) readily suggests an asymptotic (n — > oo) evaluation by the 
classic saddle-point method using |/(er)| as height function. First, the integration 
contour C is to be moved in the direction of decreasing |/(er)| as far down as possible: it 
will thus pass through some saddle-points ex* of |/(<r)|. Then for large n, 1(a) peaks near 
each of these points a*, where it makes a contribution of the order of magnitude to 
the integral: thus the highest saddle-point (s) give(s) the dominant behavior. Consistent 
asymptotic approximations can also be made inside 1(a) throughout: e.g. here, Stirling 
formulae used for T(n + const). 

This approach, for an integrand not controlled in fully closed form, partly retains an 
experimental character. We currently advocate it for this problem as a heuristic, rather 
than rigorous, tool: it predicts the global structure of the results at once, and it treats all 
the cases readily and correctly, as other techniques confirm. 

In the present problem, for large n the landscape of the function |/(<r)| is dominantly 
controlled: by its T factors, asymptotically ~ it [sin na T(2a + l)}- 1 n 2(J - 1 for finite a; and 
by the polar parts of 1(a) near its poles. The induced contour deformation starts as a 
dilation of C away from the segment [1, n] in all directions, and goes to infinity in the 
directions | arg<r| < ~ — 5. The encountered saddle-points can be of two types here (once 
n is large enough). 

1) For a on the segment (|, 1), \I(a)\ ~ n [sin na T(2a + l)]~ 1 n 2a ~ 1 Z(a) always has 
one real minimum a r (n) (tending to | as n — > oo), which will be reached by the moving 
contour; other real saddle-points lie below a — | and will not get reached here. 

2) Complex saddle-points may enter as well, for which we may focus on the upper half- 
plane alone: the lower half-plane will give complex-conjugate ("c.c") contributions. As 
long as the moving contour stays inside a half-plane {Re a > | +e}, the integrand can be 
decomposed as / = J2 k according to (HJ); then for each individual term and within the 
Stirling approximation for the T-ratio, the saddle-point equation is = ^ log |ifc(cr)| ~ 
log(a 2 — n 2 ) — 21og2a — logXfc, yielding the saddle-point location 

a k (n) =n\jlT k . (13) 

Thus any zero on the critical axis (r k real) yields a purely imaginary a k (n), not eligible: 
it lies outside the domain of validity of (0J), and its contribution would be subdominant 
anyway. So in the end, this paragraph excludes the real r k . 

The discussion then fundamentally splits depending on the presence or absence of zeros 
off the critical axis. 

[RH false] 

If there is any zero (| ± ir k ) off the critical axis, we select argr^ > and assume the 
case of a simple zero for argument's sake. Paragraph 2) above fully applies to each such 
zero: the complex saddle-point a k (n) given by (|13|) lies inside the domain of convergence 
{Re a > |} as soon asn > | Im l/r^ 1 , and for n — > +oo it gives an additive contribution 
~ \(r k + i/2) / (r k — i/2)] n (in the usual quadratic approximation of \ogI(a) around a k (n)), 
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which grows exponentially in modulus and fluctuates in phase; it will indeed exponentially 
dominate the contribution of the real saddle-point a T (n), to be computed later. 

This result can also be confirmed rigorously and more directly: by a conformal mapping 
0, the function f(z) in (JJJ) has precisely the points Zk = (r^ — i/2){jk + i/2) _1 and z% as 
simple poles of residue 1 in the unit disk; then a general Darboux theorem [6, chap. VII 
§2] applies here to the poles with \zy\ < 1, implying that the Taylor coefficients of / 
(namely, the A n ) indeed have the asymptotic form 

A ™ ~ Yl z k n + c - c - ( mod °( e£ ") > 0), n -> oo , (14) 

{|*kl<i} 

which now holds for multiple zeros as well, counted with their multiplicities. Concretely, 
A n then oscillates between exponentially growing values of both signs. 

Infinitely many zeros (| ± ir^) off the critical axis are perfectly admissible here: their 
Zk satisfy \zk\ < 1, Zk — > 1, and the corresponding infinite sums J2k z k~ n s ^ m define valid 
n — > +oo asymptotic expansions. On the other hand, the general Darboux formula (|14|) 
hopelessly breaks down if the infinitely many z^ n have identical and dominant modulus, 
which is precisely realized in the case [RH true], with all the z^ on the unit circle! 

[RH true] 

Here, Darboux's theorem only tells that A n = o(e en ) Vfe > 0, but it fails to give any clue 
as to an explicit asymptotic equivalent for A n . By contrast, the saddle-point treatment 
of the integral (|12|) itself remains thoroughly applicable. Simply now, all the are real, 
Z(a) = 0(Z(Re a) | Im cr| -3 / 2 ) in {Re a > |}, and the contour C can be freely moved 
towards the boundary {Re o = -} without meeting any of the crfc(n) (all of which are 
purely imaginary). Hence the only dominant saddle-point is now cr r (n) G (|,1); it is 
shaped by the double pole of Z(cr) at \ (itself generated by the totality of Riemann 
zeros), so that o~ r (n) ~ | + j^^. This saddle-point is however non-isolated (it tends to 
the pole), so the standard saddle-point evaluation using the quadratic approximation of 
log/(er) around <7 r (n) works very poorly. Here, it is at once simpler and more accurate 
to keep on deforming a portion of the contour C nearest to a — \ until it fully encircles 
this pole (now clockwise), and to note that the ensuing modifications to the integral are 
asymptotically smaller. Hence for [RH true], A n is given (mod o(n)) by [28j 

A n ~ (-l) n 2n Res a=1/2 1(a) = 2vn [2i2_ 2 (^(| + n) - 1 + 7) + R- X \ 

= 27m[2 J R_ 2 (logn- 1 + 7) + R^) +0(l/n) (15) 

(with ip = r'/r estimated by the Stirling formula, 7 = Euler's constant, and using the 
polar structure (0) for Z(a)). 

Prior to using (jUJ) to fix the R-j, the argument covers zeros p of a more general 
(arithmetic) Dirichlet series L(s): as long as the latter has a meromorphic structure and 
functional equation similar enough to ((s), its secondary zeta function Z{a) keeps a double 
pole at a = | j^ZJ- A related but more concrete requirement can be put on the function 
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N(T), the number of zeros of L(s) with < Im p < T: we ask that for some constants 
R-2, R-i and some a < 1 (all now depending on the chosen L-series), 

N(T) = 2T [2_R_ 2 (logT — 1) + R-i] + 5N(T), 5N(T) = 0(T a ) for T -> +oo (16) 

(implying R_ 2 > 0). If both conditions © and ([16)1 hold, then the polar coefficients of 
Z(a) in (jSJ have to be precisely the R-j from pfij) . All of that is realized in many cases 
including, but not limited to, Dedekind zeta functions ^Sl EH HZ] and some Dirichlet 
L-functions [HI 1211 EE]; see also [Hj (discussed at end); in all those instances, 5N(T) = 
O(logT). For the corresponding A n , our saddle-point evaluation then always yields: either 
(fTH|l =>- fTK|) if all the zeros have Re p = | - or the immediately general result ((HJ) 
otherwise. 

Like (fHj) before, (fT5|) can be derived quite rigorously but by still another method, 
previously unknown to us, and written for the Riemann zeros by J. Oesterle (pri- 
vate communication). We thank him for allowing us to repeat his argument here; we 
actually word it in the more general present setting (and slightly streamline it). When 
all the zeros lie on the critical line, first transform the summation (Q) into a Stieltjes 
integral (where 6{T) = 2 arctan(l/ 2T)): X n = 2 / °°[1 - cos n8(T)} diV(T), then inte- 
grate by parts: n -1 A n = 2 C sin 72$ iV(~ cot §) d6>. Now replace N(T) by its large-T 
form ()16|) neglecting SN(T) and other 0(8~ a ) terms: the error is o(l) by the Riemann- 
Lebesgue lemma, mainly because Mb}) makes SN(^ cot |) integrable over the closed interval 
[0,7r]; hence n _1 A n = J^sinn^ |[8i?_ 2 (log| — 1) + 4iiLi] d6> + o(l). Next, change vari- 
able: n9 = t, then replace upper t-bound mr by +oo again with an o(l) error; thus, 
A n = n J °° s 2i[8i?_ 2 (log j — 1) + 4i?_x] dt (mod o(n)), and the last integral evaluates 
in closed form [§] to yield (JT3J). Unfortunately, we do not see how to extend this purely 
real-analytic argument to include the [RH false] case as well. 



Recapitulation 

As we ended up with two mutually exclusive large- n behaviors for the \ n , (fTlj) and (JT3J), 
together they provide a sharp equivalence result. For the Riemann zeros, using the explicit 
values ©: 

Theorem (asymptotic criterion for the Riemann Hypothesis). For n — > +oo, the 
sequence A„ built over the Riemann zeros follows one of these asymptotic behaviors: 

- [RH true] <^ tempered growth to +oo, as 

A n ~ ~n(logn — 1 + 7 — log27r) (mod o(n)); (17) 

- [RH false] <^ non-tempered oscillations, as 

A " ~ E {v^Ti}^ - - Hod o(e £ ") Ve > 0). (18) 

{argr fc >0} Tfc ^ 

This comprehensive asymptotic statement [2H] is new on the [RH false] side to our 
knowledge, and it also completes some earlier results in the [RH true] case. 
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Our main end formula ((JUj), assuming RH) had actually been displayed by Keiper [T21 
Equation (37)], but under a somewhat misleading context: in his words, (fT7|) required not 
just RH but also "very evenly distributed" zeros, and was "much stronger than" RH; no 
details or proofs were ever supplied. All that hardly points toward our present conclusion 
that (fTTj) and RH are strictly equivalent. 

Oesterle had a proof of the statement [RH true] =>- (fTTj) (see above), but he neither 
published nor even posted his typescript [23] . 

In [21 Cor. 1(c)], rather weak exponential lower bounds A n > — ce en were shown to 
imply RH; the backward assertion [RH true] •<= (|T7jl is thus also implied by [2] (but 
cannot be inferred therefrom, as [2| never alludes to asymptotics regarding the A n ). 

Our saddle-point approach also handles both cases (JT7j) - (|18J) at once for the first time. 

Numerical data ^21 120] agree well with (|17|) for n < 7000 (and even better in the mean 
if we add the contribution like (|15|) but from the next pole of /(cr), 

S\ n = {-l) n 2n Re S(T= o 1(a) = 2Z(0) = +7/4 (19) 

[23 Equation (41)], although this correction should not count asymptotically, dominated 
as it seems by oscillatory terms). Yet the above numerical agreement is inconclusive 
regarding the Riemann Hypothesis: any currently possible violation of RH would yield a 
deviation from ()17j) detectable only at much higher n (see end of Sect. 4). 



4 An even more sensitive sequence. 

A slightly stronger difference of behavior follows for the special linear combinations ()20|) 
below of the coefficients r) n themselves (defined by (jTUj) above). 

Indeed, the definition H(s) = T(s/2)n~ s ^ 2 s(s — l)C(s) substituted into ((7j) readily yields 
a decomposition A n = S n + S n , where [21 thm 2] 

S " = ~ £ (") ( 2 °) 
is the contribution of (s — 1)C( S ) 5 an d 

S n = l- l ° gAn + 1 n + S n , with S n = ^( n \-iy(l-2-i)C(j), (21) 

is the contribution of the remaining (more explicit) factor. The large-n behavior of the 
sum S n is also computable, by the same route we followed from to (fTSj) through (fT^j) . 
First, 

integrated around the poles a = 2, . . . , n of J(cr); hence for n — > +oo, S n is asymptotically 
(-l)"- 1 ^! Res a=1 J (a) = \nty){n) + log 2 - 1 + 2 7 ], (23) 
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now (mod o(n N ) WN > 0) because J (a) has no further singularities; so that finally, using 
the Stirling expansion for %p{n) in terms of the Bernoulli numbers B 2 k, 

— 3 °° B 

S n ~ |n(logn — 1 + 7 — log27r) + - — "^k^"^ ^ n ~^ 00 ^ unconditionally. (24) 

k=l 

This at once confirms two empirical conjectures made by Maslanka 20, Equations (2.5), 
(2.8)]: mod o(n), the sequence {S n } expresses the "trend" ()17j) obeyed by the sequence 
{A n } under [RH true]; and to all orders in n, {S n } has the asymptotic expansion (|24p . 
But whether RH holds or not it makes sense to withdraw the fixed {S^l-contribution 
from the previous formulae (fl~""j) and ("T"""|) . as we did in [2E] to find: 

S n = o(n) [RH true] (25) 

further discussed in j2H I2H Hj ) , versus 



S n ~ K ~ V f Tfc + l, V + cc. (mod o(e £n ) We > 0) [RH false], (26) 

^ — 4 V Tfc — i / 2 / 

{argr fe >0} ft ' 



which gives oscillations that blow up exponentially with n. 

Still, in absolute size, any contribution from (|2fi|) will stay considerably smaller than 
(|2*3|l (the background from the set of real Tj) up to n w min{ argTfc> o}{| Im l/rfc| -1 }: i.e., 5* n 
can only reliably signal a zero violating RH up to a height | Im p\ < \Jn/2; so vice- versa, 
since such zeros are now known to require | Im p\ > 10 9 , they could only be detected by 
S n for n > 10 18 (see also [U HH], [H p. 5]). 

Our asymptotic criteria may thus not surpass others in practical sensitivity, but their 
sharpness is theoretically interesting. For instance, Li's criterion A n > (Vra) is now 
strengthened in that, beyond a finite n-range, it refers to ever larger amplitudes either 
way, i.e., delicate borderline situations such as A n — > are ruled out. 



5 More general cases: a summary 

For the more general zeta-type functions satisfying (jHJ) and pfij) . the asymptotic alterna- 
tive for the associated A n is: either (JT5j) if the generalized Riemann Hypothesis (GRH) 
holds, or (|TH|) otherwise. In the former case, the connection ()16j) =^ (JT5|l often makes the 
asymptotic form of A„ fully explicit without any calculation, and it also ensures _R__ 2 > 
in line with the corresponding generalized Li's criterion [2*1 FLU ITS! I19j. 

The A n have also been generalized to L-functions defined by Hecke operators for the 
congruence subgroup T (N) [H][IH1 specially remark 5.4]. 

More recently and in a broader setting (the A n for automorphic L-functions), La- 
garias presented an alternative approach to estimate the A n with greater accuracy, mod 
0{\fn\ogn) under GRH [T3j: in the notations of (f*?0jl -([2*T jl . he directly proves that S n 
obeys (fTHj) mod 0(1) unconditionally (thm 5.1), then that S n = 0(y/n\ogn) under GRH 
(thm 6.1). We note that his leading A n -behavior remains tied to the large-T behavior of a 
counting function according to the rule (fT6*|) ^> ffTB*)) [T4*[ Equations (2. 11-12), (1.18), (5. 2)]. 
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